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ABSTRACT 

I review my new method for solving general 1-matrix models by ex- 
panding in iV -1 without taking a physical continuum limit. Using my 
method, each coefficient of the free energy in the genus expansion is ex- 
actly computable. One can include physical information in a function 
which is uniquely specified by the action of the model. My method gives 
completely the same result with the usual one if the physical fine tuning 
is done and the leading singular terms are extracted. I also calculate in 
the genus three case and confirm the validity of my method. 
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In this review article, I want to explain a new method for solving 1-matrix 
models that I have found recently Matrix models have been vigorously studied 
as an exactly solvable toy model for two dimensional gravity. Even though 1-matrix 
models are the easiest of all matrix models, only the planar limit have been 
argued mainly until the appearance of the double scaling limit ||, £|. This limit 
enables us to convert a recursive equation derived by the orthogonal polynomial 
method and large- N expansion (N is the size of matrix) to a non-linear differential 
equation for the specific heat. The solution of this differential equation contains non- 
perturbative effects which are absent in the perturbative expansion by the Feynman 
diagrams. My method gives another way to solve the recursive equation. Instead of 
taking the double scaling limit, I expand the recursive equation in iV" 1 and solve it 
order by order. What an intuitive way! However, the equations for the expansion 
coefficients of r e (r e will be explained shortly later) derived in this way is hard to 
solve even when the expansion order is low. The obstacle to advance was the fact 
that these equations have summations over paths ||. To overcome this difficulty, I 
found a beautiful relation between these summations and a function W(r) which is 
specified uniquely by the action of the model. Using this relation I find the expansion 
coefficients of free energy until genus three In appendix I will show the explicit 
result for the genus three case. In principle one can continue this calculation to 
higher genus. However, the calculation will become rapidly complicated for large 
genus. 

I follow the notations of my previous paper |J. Since in this short review I cannot 
explain all of the notations, I will sometimes omit explanation for some undefined 
signs. Consult [Q for further information. A 1-matrix model is defined by an integral 
over an N x N hermitian matrix $, 

Z N ({\}) = b$e- NS ^, 

1 oo 

5[$,{A}] = tr $ 2 + ^A p $ 2 ^ 
1 P =i 

where {A} denotes the set of coupling parameters. The main purpose of the game 
is to find the expansion coefficients of the free energy 

E N ({\}) = ^og^^ = ±N-^e^(a\{\}). (2) 

In (H) I indicated a 2 -dependence of e^ h \ This parameter a 2 is a real solution of 
W(r) = 10, where W(r) is the function I mentioned before 



W{r) 



f (2P + 2)! 
^p!(p+l)! P 



(3) 



1 If such a solution does not exist, the matrix integration will be complex and loses its physical 
meaning. 
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Using the orthogonal polynomial and the Euler-Maclaurin formula the contribu- 
tion from genus h is described as follows 



e (h) 



[ dx(l - x)f 2h (x) + K 2h 
Jo 



-E(-i)*7^{(i-*)^.)(*)} 



(4) 



(2a)! 

where B s are the Bernoulli numbers. This formula is my starting point. The uni- 
versal contribution which is survived in the continuum limit lies in the first term on 
the right hand side. In f 2s are expansion coefficients of the function 

r (r) 00 

f e (:r)=fog^ = £e*f !ta (x), (5) 

X s=0 

where e = N^ 1 . Since the boundary term K 2 h can be derived from the function 
W(r), what we need to compute is to solve a recursive equation for r € (x) 



x = r e (x) 



l + ]T2(p + l)A p £ t[r e (x + Si (l)e) 



(6) 



from which one can find the moments f 2s in (£|). In (^j) L p is a set of all paths with 
length (2p + 1) discussed in my previous paper and — p < Si(l) < p are p integers 
assigned to each path I. These so-called "height" variables are connected to the 
set of all combinations which pick up p elements out of a set of (2p + 1) integers 
{1, 2, . . . , 2p + 1}. If one of the combination is {ax, . . . , cr p } , (a% < ■ ■ ■ < a p ), each 
integer <jj is related to a height variable for a path / 6 L p according to the following 
simple rule: 

Si = Gi -2i, i = 1,2, . . .,p. (7) 

If the action of the model is suitably fine tuned and the double scaling limit is 
taken, this recursive equation turns to the string equation. Alternatively I solve (^j) 
by expanding as r e (x) = I] s e 2s r 2s (x). One should note that the equation for r is 
x = W(tq). Thus the boundary x — 1 in (Q) corresponds to tq = a 2 . This boundary 
contributes to the non-analyticity of and the another boundary x = to the 
regular part in terms of coupling parameters. 

It is convenient to introduce a number J p (m,...,n m ) specified by a Young's 
diagram (ni, . . . , n m ) whose number of rows is m. This Young's diagram corresponds 
to a partition of a positive even integer 2n, i.e. , J2 n i = 2n. These numbers play an 
important role in my method, for they appear in the expanded recursive equations 
derived from @. If f(x) is an arbitrary infinitely different iable function, the number 
is defined by a following Taylor expansion 

p 00 
£ II / + Si(t)e) = e 2n F 2n (p) , 

ieLp i=l n=0 
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F 2n (p) = J pK • • • > n m )f^\x) ■ ■ ■ f M (x){f(x)y- m , (8) 

{ni,...,n m )&2n 

where Y 2n is the set of partitions of 2n. The number vanishes when p < m. These 
numbers depend on p which parameterizes the path length. I have decided the in- 
dependence of these numbers guided by two kinds of arithmetic generating functions: 

oo 

2p+iC p (x) = H Ip(ni, ■ ■ ■ ,n m )x 2n , (9) 

n=0 (ni,...,n m )eY2„ 

oo 

T p (x) = £/ p (2n)a; 2 ". (10) 

n=0 

Remarkably these generating functions are both described by g-deformed integers, 
where q = exp(x/2). The answer for this problem can be presented in a simple way. 
To do this it is convenient to define a generating function for these numbers 

oo 

J(r ;K...,n m ))= £ 2(p + l)\ p I p (n u . . . , n m )r p - m . (11) 

The p-dependence of I p {n\, . . . , n m ) is described by the following remarkable relation 
between the generating function and the function W(tq), 

n-l 

J(r ; (m, . . . , n m ) ) = ^ ^(m, . . . , n m y o W^ +m+ ^ (r ) , (12) 

where ai(ni, . . . , n m ) are positive rational numbers. These n rational numbers are 
completely fixed by n explicit values of I p (ni, . . . , n m ) for small p. Using a computer 
these values are computable. One meets the generating functions defined in ( |TT|) 
and their derivatives with respective to r when the large-iV expansion of @ is 
performed. For example, one obtains a equation for r 2 (rather for f 2 = r 2 /r ) from 
the e 2 order, 

r 2 W(r ) + r {/(r ; (2) )r (2) + J(r ; (1, 1) )(r (1) ) 2 } = . (13) 
The relation ( TT2] ) is easily verified for the Young's diagrams (2), and (1,1): 

7(r Q ; (2) ) = \w^\r Q ) , J(r ; (1, 1) ) = ^( 3 )(r ) . 

Since one can write any derivatives in terms of derivatives of W(r ) using the 
relation x = W(r ), the moment f 2 can be described only in terms of W(ro). 

These expanded recursive equations give the solutions for the expansion coeffi- 
cients of r e . These coefficients are always written in the following form for h > 1, 
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The above form has an important meaning: the integral in (|J) can be explicitly 
performed and only the boundaries x = 0, 1 contribute to the free energy. I must 
introduce a subset X k C Y 2k , where Y 2 k is the set of partitions of 2k to explain the 
functions Q^ 1 ' in ([14]). This subset is defined as a set of the partitions of 2k whose 
number of rows of the corresponding Young's diagram are k, 

X k = {{n 1 ,...,n k )\J2n i = 2k}. (15) 

i 

Then I define sets of rational polynomial functions. If f(x) is an infinitely differen- 
tiable function, the sets are defined as follows 



<?*,.(/(*))= Q{f[f ini) }/(f) 

(m,...,n k )eX k U=l J 



(16) 



Then the functions Q^(r ) are elements of such sets, 

nP(r ) E Q 2{h - 1)+sMh - 1)+s (W(r )) • (17) 

The actual computation process of e^ h > is quite boring, so I show only the result. 
Since the spherical contribution is well-known, I omit this case. The genus one 
contribution is surprisingly simple and is exceptional comparing to higher genus 
cases, 

e U = -±log{a 2 W'(a 2 )}. (18) 

For h > 2 the explicit form of e^ 1 ' is very complicated, but they are always written 
in a following form, 

e<*>= £ Ei h \a*y* + (-l) h —?±— (19) 

a=-2(ft-l) 4A H"' l > 

where the coefficients E^ h \ro) are the same kind of rational polynomial functions of 
derivatives of W(r ) as fi^ except that the subscript s might be negative, 

Ei h) (r ) e Q 2 

(h-l)+s,4(h-l)+s 

(W(r )) . (20) 

Here the coefficient with the greatest s is equivalent to fij^i- The Contributions 
from negative s are understood through a relation, 

j ^vy» = (-i)'^P (^P - 2 - (21) 

Furthermore Eq 1 ^ can be constructed from the following relation, 



4 h) (0) = -hm 

r— >0 



E S?)(r)r- + (-r fc L ' 



s=-2(h-l) 



4/i(/i - 1) H/(r) 2/l " 2 



(22) 
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The most important thing to note in (|T9|) is that for h > 1, the x = boundary 
contributions are offset and there is no term that depends on the derivatives of 
W(r ) at r = 0. The origin of this cancellation of boundary contributions is still 
mysterious. In appendix I concretely compute e^ 3 -*. The reader will understand how 
complicated the actual form of the free energy is. 

One must fine tune the action to treat the 1-matrix model defined by (fj) as a 
model of two dimensional gravity. The critical exponent and the physical observables 
are included in this model by the fine tuning. I think that the most heuristic way 
for the fine tuning is simply described as follows 

m 

Swp^r}] = N + ^ Tifffim . (23) 

i=-2 

In ( |2"B| ) Pn-kW\ expresses a discretized scaling operator. The i = term corre- 
sponds to the cosmological term, (n ,M ,r ) = (0, 0, t) . The terms with negative i 
specify the criticality of the system, (n_i, M_i, r_i) = (k, —1,6), (n_2, M_ 2 , r_ 2 ) = 
(k, — 2,— 1). The integer k is related to the string susceptibility 7 = —1/k in the 
continuum theory. The terms with positive i describe source terms for scaling op- 
erators and rii/k is equivalent to the gravitational dimension of the i-th scaling 
operator. The superscripts Mi and b are non-universal counterparts of rii and k 
respectively. Since the formal gravitational dimension for the (— l)-th scaling op- 
erator is 1, this term corresponds to a marginal operator and b can be interpreted 
as a marginal deformation parameter. Though I do not show the explicit form of 
the discretized version of the scaling operators, it is sufficient to show the converted 
function WkfiiT, { r }) derived from the fine tuned action (|23|). 

«'*W) = i-Ey (i-jqpj) n. (24) 

Using this function, one can control even the non-universal information systemati- 
cally. For example, one can find that the regular part of the spherical free energy is 
at most a bilinear function with respective to the couplings 73, (i = 0, 1, . . .). 

Now I fine tune the action to handle the free energy (Q) as a generating function 
for correlation functions of physical observables. From the physical point of view a 
continuum limit must be taken. The rule for extracting the universal terms from 
the result of ( |T^ ) is quite simple: replace the derivatives of W(r ) at r = a 2 in 
the term £ , ^ 1 (a 2 )a 2< - /l_1 - ) according to the following manner, 

where fo(t) is the solution of t — /q — YA^fo*' Note that at this stage the non- 
universal parameters Mj and b disappear. When all sources are turned off, this 
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function fo(t) corresponds to the specific heat in the continuum limit. Finally the 
universal contribution e h \t) becomes a rational function of the derivatives of fo(t), 

e { o h \t) e Q 3M) 4 M (./o(t)) , (h>2). (26) 

One can compute any correlation functions of physical observables from this uni- 
versal function e^\t). Their t-dependence completely agrees with the prediction of 
the continuum theory. When all source terms are absent, one can confirm that the 
t-dependence of the free energy actually coincides with the result of the Liouville 
theory 0, 

e^(t)| M=o oct 2 -^), (27) 

where, 

2-Mh) = i.3(/i-l)-6(/i-l)-4(/i-l)(i-l 

= (1-/0(2-7). (28) 

In this way my method reproduces all of the known results. 

At the end of this short article, one should note that the matrix size N is still 
finite. Discarding except for the most singular terms corresponds to taking the 
double scaling limit, N — * oo, t — > and NH 2-1 — > const. 



Appendix 



In this appendix I give detailed results for genus three. The starting point for genus 
three is read from (HI), 



/ dxil — x)r 6 + Kq 
Jo 

+_L ((i _ x)h f) 1 _ _L ((i _ x)h f) 



o 720 

+ —^—((l-xy o f ] \ (29) 
30240 VV ; ; o ' v ; 

where K e is given by derivatives of W(r ) at r = 0, 

K 6 = — ^[3714iy (2) (0) 5 - 6340iy (2) (0)V (3) (0) + 1823iy (2) (0)iy (3) (0) 2 
5760 

+i455iy (2) (o)V (4) (o) - 3i8iy (3) (o)iy (4) (o) - 2iow (2) (o)iy (5) (o) 

+151^ (6) (0)]. (30) 

The solutions for f 2 = r 2 /r and f 4 = r 4 /r — (^2/ r o) 2 /2 were derived in my previous 
work. Thus the only unknown function in (p9| ) is 

r 6 r 2 r 4 1 /r 2 \ 3 
r 6 = — - — — + a — • ( 31 ) 



ro r r 3 \r 
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To find this function according to my procedure we need the generating functions 
for the partitions of until six defined in (|ll|). One should rely on a computer to do 
this. Since the generating functions for the partitions of two and four have been 
implicitly derived in my previous work, I show here only the functions for eleven 
partitions of six 2 : 



7(6 
7(5,1 
1(4,2 
1(3,3 
7(4,1,1 
7(3,2,1 
7(2,2,2 
7(3,1,1,1 
7(2,2,1,1 
7(2,1,1,1,1 



7(1,1,1,1,1,1 



2160 



-W {2) 



1 ^(3) + -L r lw^ 



360 



r l 



-T , 

840 



720 840 280 



23 

1440 
17 

2160 
31 



W {3) + 



17 

560 
11 



r 



19 

2520 ? 
23 



W {3) + —r W {4) + ^r^r 2 n W {b) 



630 



6048 



5040 
5 



1008 



W {6) 



1008 



1260 7560 
11 t j7- f 4) 227 
560 10080 



W (5) + 



2520 
61 

15120' 



24 432 252 

J^(5) + ^ ro iy(6) + 
288 576 10080 



rniy (7) 



32 



181 
10080 



r 



L^(7) + J_ roH f(8) + 



448 



960 



17 
8640 ? 



10368 



(32) 
(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



Then it is possible to solve the recursive equation for r 6 and the results are indeed 
written in the form of 



1 



d 



W'(r)dr Q J [ 



r Q\' } + r 2 Q ( 2 



(43) 



Explicit form of the coefficients Vlf^ obeys the rule (117]) 



(3) 



[18039(l^ (2) ) 4 - 2A312W'(W (2) ) 2 W {3) 
+3182(PT) 2 (W /(3) ) 2 + 5088(W') 2 W {2) W {4) 

-528(W /, ) 3 W /{5) ]/(725760(W^ / ) 8 ) > 

[2iooo(w /(2) ) 5 - sssssphw^) 3 ^ 3 ) 

+ 12114(^') 2 W /(2) (W /(3) ) 2 + 9A95(W') 2 (W (2) ) 2 W^ 
-2315(W") 3 W (3) W (4) - 1494(W /, ) 3 W /(2) W /(5) 



(44) 



2 Here I{n\, . . . , n m ) is an abbreviation for J(ro; (rii, . . . , n m )). 



8 



+ 119(W') i W {6) ]/(18U40(W') 9 ) , (45) 
fif = [34300(W (2) ) 6 - 810Q0W'(W {2) ) 4 W {3) 

+43050(iy') 2 (^ (2) ) 2 (^ (3) ) 2 - 2915(W) 3 (^ (3) ) 3 
+222Q0(W') 2 (W {2) ) 3 W {4) - 13A52(W') 3 W {2) W {3) W {4) 
+Q07(W') 4 (W^) 2 - A28A(W') 3 (W {2) ) 2 W {5) 
+100Q(W') 4 W i3) W i5) + 539 {W'fW^W^ 
-35{W') 5 W i7) }/(362880{W') 10 ) . (46) 

As the function f 6 is explicitly obtained, it is easy to compute the generating 
function fl29|) . It is certainly written in the form of (|19"D , 



e 



W= E^H^-ttL- ( 47 ) 



1008 

The coefficients of (|47D are very complicated but they indeed satisfy (|20|). 

Bf ) = nf ) (a 2 ). (4? 
Sj 3) = -[21420^ (2) (a 2 ) 5 -40110^'(a 2 )^( 2 )(a 2 )V (3) (a 2 ) 

+ 12783W'(a 2 ) 2 W i2 \a 2 )W (3 \a 2 ) 2 + 10170W'(a 2 ) 2 W {2) (a 2 ) 2 W {4) (a 2 ) 



2A88W\a 2 fw {3 \a 2 )W^\a 2 ) - 16MW'{a 2 fw i2) {a 2 )W {5) { 



a 



+133^'(a 2 ) V (6) (a 2 )]/(3628801^'(a 2 ) 9 ) , (49) 
E^ = [1575iy( 2 )(a 2 ) 4 - 1800iy'(a 2 )^ (2) (a 2 )V( 3 )(a 2 ) 

+200W'(a 2 ) 2 W^ 3 \a 2 f + 388W\a 2 fw {2 \a 2 )W^\a 2 ) 

-2AW'(a 2 fw {5 \a 2 )}/(725760W'(a 2 f) , (50) 
E {3 1 = [15W^ 2 \a 2 f - 10W'(a 2 )W^(a 2 )W^(a 2 ) 

+W'(a 2 fw^(a 2 )}/(Q0A8W'(a 2 ) 7 ) , (51) 

E {3) 2 = [15^ (2) (a 2 ) 2 -4^'(a 2 )W /(3) (a 2 )]/(6048^'(a 2 ) 6 ), (52) 

E {3) 3 = iy (2) (a 2 )/(504iy'(a 2 ) 5 ), (53) 

E {3 \ = l/(1008W(a 2 ) 4 ) . (54) 

The validity of this result is checked by applying the simplest function W(r ) = 
r + 12Air 2 , (Ai > —1/48). The apparent form of e^ 3 ^(a 2 ) certainly has the predicted 
form in ||. If the physical fine tuning is performed and the universal terms are 
extracted according to the rule fl2"5|), one can also confirm the validity of my results 
by comparing it to the known results derived by the usual method. 
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